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Summary 

An experimental influence coefficient technique was used to obtain unsteady aerodynamic influence 
coefficients and, consequently, unsteady pressures for a cascade of symmetric airfoils oscillating in pitch 
about midchord. Stagger angles of 0° and 10° were investigated for a cascade with a gap-to-chord ratio of 
0.417 operating at an axial Mach number of 1.9, resulting in a supersonic leading-edge locus. Reduced 
frequencies ranged from 0.056 to 0.2. The influence coefficients obtained determine the unsteady 
pressures for any interblade phase angle. The unsteady pressures were compared with those predicted by 
several algorithms for interblade phase angles of 0° and 180°. 


1.0 Introduction 

This report represents a research effort carried out nearly a decade ago when the NASA Lewis 
Research Center (now the NASA Glenn Research Center) was developing propulsion technology in 
support of the supersonic throughflow fan (SSTF) engine. A detailed description of this engine and its 
benefits, as well as associated research, is given in references 1 and 2. A unique characteristic of the 
SSTF is that the rotor operates in supersonic axial flow and hence all rotor blade bow shocks exist 
downstream of the locus of blade leading edges. A rotor with this flow configuration is characterized as 
having a supersonic leading-edge locus (SLEL). Experimental research on the SSTF concept prior to the 
work performed at NASA Lewis was extremely limited (refs. 3 to 5). Therefore, a research effort was 
begun to evaluate the concept and potential of the SSTF engine (refs. 6 and 7). During the initial design 
of the SSTF rotor, aeroelastic stability became a concern. Consequently, a computer algorithm (ref. 8) 
was developed to determine the aeroelastic stability of the SSTF. This algorithm is based on Lane’s 
(ref. 9) unsteady aerodynamic theory for an oscillating cascade of flat plates with an SLEL. The 
algorithm was incoiporated into an existing aeroelastic code (ref. 1 0) for use in determining the 
aeroelastic stability of the SSTF rotor (ref. 11). 

In an effort to improve on the aeroelastic analysis capabilities while still providing for efficient 
analysis turnaround time, the algorithm based on Lane’s theory was modified to include the effects of 
thickness and camber in the noninterference regions of a cascade with an SLEL (refs. 12 and 13). In 
addition, computational fluid dynamic (CFD) codes were being developed to predict the unsteady lift and 
moment coefficients necessary for determining aeroelastic stability (refs. 14 to 18). The desire to 
correlate the codes with experimental data resulted in an experimental effort to measure unsteady 
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pressures on an oscillating cascade with an SLEL. This effort was first described in references 19 and 20. 
This report compares the experimental data with values predicted by the algorithms of references 8,12, 
and 18. 

Experimental unsteady pressure data are presented for a linear cascade of airfoils with an SLEL 
oscillating (with amplitudes of 1.2° about 0° angle of attack) in pitch about midchord having a gap-to- 
chord ratio of 0.417 and stagger angles of 0° and 10°. Experimental influence coefficients were obtained 
that determine unsteady pressures for any interblade phase angle. Specifically, unsteady pressure 
coefficients for interblade phase angles of 0° and 180° were determined from the experimental influence 
coefficients. 

The influence coefficient technique used to acquire the unsteady pressure data is described next, and 
then the experimental facility and the data acquisition and reduction techniques. Experimental results for 
stagger angles of 0° and 10° and for interblade phase angles of 0° and 180° are then presented. Although a 
stagger angle of 0° would not typically simulate an actual turbomachinery situation, it nonetheless 
simplifies the experimental process (as made clear subsequently) and at the same time provides 
experimental data for correlation purposes. 


2.0 Influence Coefficient Technique 

2.1 Preliminaries 

A linear influence coefficient technique based on the work of reference 2 1 was used to compute the 
unsteady pressure distribution on the oscillating cascade. This technique determines unsteady pressure 
coefficients C p for any interblade phase angle p. The complex-valued unsteady pressure coefficient is 
defined as 


r = 

^ P 


9oo^l 


( 1 ) 


where P\ is the complex -valued first harmonic of the airfoil surface static pressure, q fj is the free-stream 
dynamic pressure, and ai is the amplitude of airfoil oscillation (radians). All symbols are defined in 
appendix A. 

For a finite cascade of 2 N + 1 airfoils executing constant harmonic oscillations with a constant 
interblade phase angle P, the unsteady pressure coefficient C p , can be expressed as a Fourier series of 
complex -valued unsteady aerodynamic influence coefficients C n (ref. 21) 

Pu,l 

N 

C P J*. P)= z c; m*"* (2) 

n=—N 

The superscript is the index number of the airfoil that is oscillating, and the subscripts u and / refer to the 
upper (suction) and lower (pressure) surfaces of the reference airfoil ( n = 0), respectively. 

Figure 1 illustrates a cascade of airfoils. In a cascade of airfoils with an SLEL, each interblade 
passage is aerodynamically isolated from all other interblade passages. Therefore, in addition to the 
reference airfoil (n = 0), only the two airfoils n = — 1, 1 that bound the reference airfoil contribute to the 
unsteady pressures on the reference airfoil. In other words, only airfoils with indices n = — 1,0, 1 
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contribute to the unsteady pressure on the reference airfoil. More specifically, both surfaces of the 
reference airfoil (n = 0), the lower surface of the upper airfoil (n = 1), and the upper surface of the lower 
airfoil ( n = -1) contribute to the unsteady pressure coefficients on the reference airfoil. Therefore, 
equation (2) becomes 


= £ c;jx)e' 


infi 


n=—\ 


(3) 


Figure 2 illustrates a cascade of airfoils with an SLEL. Portions of each airfoil surface are isolated 
from any aerodynamic influences originating from neighboring airfoils. These surfaces are denoted as 
noninterference surfaces. In particular, the unsteady pressure distribution on the noninterference surfaces 
of the reference airfoil are not influenced by the presence of airfoils n = — 1, 1. Therefore, the unsteady 
pressures on these surface regions of the reference airfoil are not influenced by interblade phase angles. 
All influence coefficients on the noninterference airfoil regions are null except for the influence 
coefficient due to the motion of the reference airfoil itself. Thus, 

C” (x,p) = 0 n± 0 x < xf l (4) 

Pu,l 

where x u s and x\ are the shock impingement locations on the upper and lower surfaces of the reference 
airfoil, respectively (see fig. 2). 

The unsteady pressure coefficients upstream of the first shock wave impingement location for the 
upper and lower reference airfoil surfaces are 

C Pu (x,fi) = C° Pu x<x“ (5) 

and 

C p , (tP) = C® x<x[ (6) 


respectively. The unsteady pressure coefficients downstream of the first shock wave impingement 
location for the upper and lower reference airfoil surfaces are 

(7) 

( 8 ) 


C (x,|3) = C° +C l ^ x>x l ‘ 

^ u Fu Fu 

and 

C pi (x,V) = C° pi+ C-y® x>x{ 


respectively. 

Equations (5) to (8) may be combined by introducing the asymmetrical unit step function defined as 
(ref. 22) 


U_(x-x s ) 


f 0 x < x s 
[l x > x s 


(9) 
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Combining equations (5) and (7) and incorporating equation (9) yield 


C„ (jc,P) = C° +C 1 e®U_(x-x l ‘) (10) 

Pu y Pu Pu V s 

Also combining equations (6) and (8) and using equation (9) give 

c Pt ( X >P) = Cp t +C~yv U_(x-x l s ) (11) 


Equations (10) and (11) can be used to determine the unsteady pressure and consequently the unsteady 
lift and moment coefficients for a cascade with an SLEL. 

The two-dimensional unsteady lift and moment coefficients are defined as 


and 



(12) 

(13) 


respectively, where c is the airfoil chord and x is the streamwise coordinate. The unsteady pressure 
difference coefficient A C p for the reference airfoil is defined as 


AC -C -C 

^ P ^ Pi ^ Pu 


(14) 


Upon substitution of equations (10) and (11) into equation (14) the unsteady pressure difference 
coefficient becomes 


A C„=C° +C~ X e~^U_(x-x l s )-C° -C l «?%_(*-*“) 


Pi 


Pi 


Pu 


Pu 


(15) 


2.2 Techniques 

Figure 3 shows five techniques or scenarios for measuring the unsteady pressure influence 
coefficients by using a cascade of two or three airfoils. The first technique is accomplished in three steps 
(fig. 3(a)). It uses three airfoils in the first step with only two airfoils required in the second and third 
steps. In each step a different airfoil is oscillated while unsteady pressures are measured on the reference 
airfoil. The reference airfoil is instrumented on both its upper and lower surfaces and remains in the 
middle location of the cascade for each step. The airfoil oscillation system for this technique must be 
able to oscillate an airfoil in each of the three different airfoil locations within the cascade. 

The second technique (fig. 3(b)) uses three airfoils in the first step, with only two airfoils required in 
the remaining steps, just as for the technique described above. The reference airfoil is instrumented on 
both its upper and lower surfaces and is located in a different position in the cascade for each step. 
However, unlike the first approach the middle airfoil is oscillated in each step while unsteady pressures 
are measured on the reference airfoil, which is in a different location for each step. The airfoil oscillation 
system for this technique need only oscillate the airfoil in the middle location of the cascade. 
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For symmetric airfoils the technique illustrated in figure 3(c) may be used and was the method 
employed in this work. This technique uses a cascade of two monoconvex airfoils instrumented on their 
convex surfaces. It requires two steps, each with equal and opposite stagger angle. This technique will 
work for nonsymmetric airfoils, where one airfoil is instrumented on its pressure (lower) surface and the 
other airfoil is instrumented on its suction (upper) surface. Nonsymmetric airfoils must be inverted in 
step 2 as indicated by the coordinate system orientation and must be switched in location as well. The 
airfoil oscillation system for this technique need only oscillate the lower airfoil of the cascade. This 
technique will be described in detail shortly, but for now two additional techniques are presented. 

Figure 3(d) illustrates an additional technique using a two-airfoil cascade. This technique consists of 
four steps and uses both positive and negative stagger. One airfoil is instrumented on both surfaces. In 
each step the lower airfoil in the cascade is oscillated. The instrumented airfoil is in a different position 
in the cascade from step to step. For symmetric airfoils this technique can use monoconvex airfoils, with 
their convex surfaces facing each other. Nonsymmetric airfoils must be inverted in steps 3 and 4 as 
indicated by the respective coordinate system orientations. The disadvantage of this approach is the need 
to perform four wind tunnel runs. 

Finally, figure 3(e) presents a technique using three airfoils. The middle airfoil in the cascade is 
instrumented on both surfaces. The upper and lower airfoils are instrumented on their lower and upper 
surfaces, respectively. All influence coefficients are obtained from the four instrumented surfaces in one 
step wherein the middle airfoil is oscillated. 

As already stated, the technique illustrated in figure 3(c) was used in this work and is denoted as the 
two-airfoil technique. This novel technique has several advantages over the others presented. The two- 
airfoil technique requires only two steps, as opposed to the three steps required for the techniques 
illustrated in figures 3(a) and (b) and the four steps required for that illustrated in figure 3(d). Although 
only one step is required for the technique illustrated in figure 3(e), the disadvantage of this approach is 
the need to instrument four surfaces and the resulting increase in data acquisition efforts to sample all 
four surfaces simultaneously. 

Although not used in this work, the three-airfoil technique illustrated in figure 3(a) is presented next 
to facilitate a better understanding of the two-airfoil technique, which will be presented subsequently. 

2.2.1 Three-airfoil technique (one instrumented airfoil). — With the upper and lower surfaces of the 
reference airfoil instrumented (see fig. 3(a)), all coefficients corresponding to a given Mach number, 
reduced frequency, and stagger angle can be obtained in three steps as follows: 

1. Oscillate the reference airfoil (n = 0) and measure the unsteady pressure on the reference airfoil 
to obtain C® and C'jj as shown in step 1 of figure 3(a). 

2. Oscillate the lower airfoil (n = -1) and measure the unsteady pressure on the lower surface of the 
reference airfoil to get c j as shown in step 2 of figure 3(a). 

3. Oscillate the upper airfoil (n = 1) and measure the unsteady pressure on the upper surface of the 
reference airfoil to get C 1 as shown in step 3 of figure 3(a). 

Pu ~ 

2.2.2 Two-airfoil technique (two instrumented airfoils). — In order to acquire the four influence 
coefficients for a cascade with stagger, two steps are necessary as shown in figure 3(c). Each step would 
normally correspond to a separate wind tunnel run: one run with the airfoils at positive stagger and the 
other with the airfoils at an equal negative stagger (fig. 3(c)). The reduced frequency must be the same 
for both wind tunnel runs. For a case with stagger all steps may be accomplished in one wind tunnel run 
if some mechanical arrangement is used to automatically adjust from positive to negative stagger during 
the same wind tunnel run. For the case of 0° stagger angle, only one wind tunnel run is necessary, as 
discussed in section 2.2.3. 
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In the two-airfoil technique two airfoils are instrumented on their convex surfaces. The lower and 
upper airfoils are designated as airfoils A and B, respectively. Influence coefficients are acquired for a 
cascade with stagger in the following manner: 

1. With the airfoils at positive stagger as shown in step 1 of figure 3(c), airfoil A is oscillated. As 

shown, influence coefficients C° and C _1 are obtained in the same wind tunnel run. Airfoil A 

Pu Pi 

simultaneously simulates the reference airfoil and the lower airfoil (n = -1) of the three-airfoil 
technique as shown in steps 1 and 2 of figure 3(a). Airfoil B simultaneously simulates the upper 
airfoil (n= 1 ) and the reference airfoil of the three-airfoil technique as shown in steps 1 and 2 of 
figure 3(a). 

2. With the airfoils at negative stagger as shown in step 2 of figure 3(c), airfoil A is oscillated 
again. Influence coefficients and (A are obtained. Airfoils A and B simulate the upper 
airfoil {n= 1 ) and the reference airfoil, respectively, of the three-airfoil technique as shown in 
step 3 of figure 3(a). In this case the pitch-up (clockwise rotation) of airfoil A in step 2 (fig. 3(c)) 
represents a pitch-down (counterclockwise rotation) of airfoil 1 in the three-airfoil technique 
(fig. 3(a)). 

2.2.3 Two-airfoil technique (zero stagger angle ). — When the two-airfoil cascade has 0° stagger 
angle, equation (15) can be further simplified. At 0° stagger angle the shock wave impingements are at 
the same x/c locations for the upper and lower surfaces of airfoils A and B (generally only true for 
symmetric airfoils or, in this case, monoconvex airfoils with their convex portions facing each other). 
Hence, 

x s = 4 = x s ( 16 ) 


A pitch-up (clockwise rotation) of airfoil A in step 1 of figure 3(c) corresponds to a pitch-down 

(counterclockwise rotation) of airfoil A in step 2 of figure 3(c). Thus, the motion of airfoil A in step 1 of 

figure 3(c) is 180° out of phase with the motion of airfoil A in step 2 of figure 3(c), from the viewpoint of 

the cascade coordinate systems shown. Because all influence coefficients are acquired with respect to the 

oscillatory motion of airfoil A, it is evident that C -1 will be 180° out of phase with respect to C 1 . 

Pi — i — i P U 

However, their magnitudes will be equal. For this case of 0° stagger, coefficients C~ and C p are 

measured on airfoil B in the same wind tunnel run. Therefore, 


C _1 = C 1 e in = -C l 

Pi Pu Pu 


(17) 


Likewise, coefficients C° and C° are measured on airfoil A in the same wind tunnel run. It can also be 

Pu Pi - 0 - 0 

seen from steps 1 and 2 of figure 3(c) that coefficients C p and C p are 180° out of phase and have the 
same magnitudes. Therefore, 


C° =C° e in = -C° 
Pu Pi Pi 


(18) 


Substituting equations (16) to (18) into equation (15) yields 

AC„ = 2 C° - C 1 ( e~ * P + e' P ) U_(x-x„) 
P Pi Pu\ ! v sJ 


(19) 
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Let C pA and C represent the influence coefficients measured at 0° stagger angle on airfoils A and B, 
respectively, resulting in 


and 


r =C° =C° e in 
Pu Pi 


c v = C _1 = C 1 e in 

P B Pi Pu 


= -c u 


Pi 


= -C L 


Pu 


( 20 ) 

( 21 ) 


Substituting equations (20) and (2 1 ) into equation ( 1 9) gives 


\C p =-2Cp A+ C pB (e-V + eV) U_(x-x s ) 


( 22 ) 


Equation (22) can be used to calculate the unsteady pressure difference coefficient for a cascade of 
symmetric airfoils with an SLEL having a 0° stagger angle. 


3.0 Experimental Facility 

3.1 Wind Tunnel 

The experimental oscillating cascade facility (OCF) used to obtain the unsteady pressures was 
conceived as a modification to the existing Ohio State University Aeronautical and Astronautical 
Research Laboratory’s high-Reynolds-number blowdown wind tunnel. The modification consisted of 
replacing the transonic tunnel section with a supersonic tunnel section. The OCF (fig. 4) consists of a 
settling chamber, a supersonic tunnel section, an exhaust pipe, and an oscillating mechanism. The 
complete oscillating mechanism is not shown in figure 4. The settling chamber is approximately 
10 ft long and 3 ft in diameter. The tunnel is 58 in. long with a constant-area test section 6 in. wide by 
1 2 in. high. 

Referring to figure 4, the dry air enters the dual-can squirter and is diffused radially as it enters the 
settling chamber. The dual-can squirter consists of two concentrically mounted cans with solid 
downstream ends and porous sides. Upon entering the settling chamber the airflow then travels through a 
porous plate and through a section of 9-in. -long honeycomb. The honeycomb straightens the airflow by 
removing radial velocity components. After leaving the honeycomb flow straightener the air passes 
through three screens which dampen any remaining turbulence in the flow. After passing through the last 
settling chamber screen the flow then enters the convergent portion of the wind tunnel, passes through 
the throat, and is expanded to supersonic flow in the test section. The flow exits the facility through a 
circular pipe and is exhausted into the atmosphere. 

Figure 5 shows a schematic diagram of the air supply system for the blowdown OCF. Two storage 

3 

tanks provide 1500 ft of dry air at 2400 psi. The air is compressed and pumped through an air dryer into 
the tanks by two four-stage compressors. A central valve system controls the airflow to the OCF. The air 
is admitted to the stagnation/settling chamber by opening a 1/4-tum plug valve. The flow rate is set by a 
cone valve. A manually operated butterfly valve connects the OCF to the high-pressure air supply. 

Figure 6 shows the aluminum wind tunnel section with the north side wall and north nozzle block 
removed. The tunnel sides consist of the nozzle blocks (fig. 7) and the side walls. The roof and floor are 
each one contiguous piece as shown in figure 6. Fasteners join the nozzle blocks and side walls to the 
roof and floor. The exhaust flange (fig. 6) connects the tunnel to the exhaust pipe. The airfoil windows 
mount in the nozzle blocks as shown in figures 6 and 7. The windows may be rotated in the nozzle 
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blocks, positioning the upper airfoil either upstream, downstream, or directly above the lower airfoil and 
thus providing various stagger angles. The window angular orientation and thus cascade stagger angle 
was determined by measuring the inclination from horizontal of a flat plate placed on two dowels 
extending outward from the window as shown in figure 8. The windows were secured in place by using 
four toe clamps also shown in figure 8. 


3.2 Airfoils 

Figure 9 shows one of the two airfoils that compose the supersonic cascade. The airfoils are 
monoconvex with an instrumentation cover machined into the outer surface (away from the cascade 
passage) as shown in figure 10. Figure 1 1 shows the airfoil with the instrumentation cover installed. 
Figure 12 shows the cascade with one nozzle block removed. The airfoils were machined from 17-4PH 
stainless steel and heat treated to improve fracture toughness. A circular trunnion at each spanwise end of 
the airfoils was followed by a rectangular boss. The circular trunnion fit into a sintered bronze bushing 
mounted in the tunnel window. This trunnion supported the aerodynamic loading while allowing angle of 
attack orientation and/or oscillatory motion of the airfoil. The boss on the south end of the lower airfoil 
(airfoil A) provided the attachment point for the oscillating mechanism, while the boss on the north end 
provided a reference surface for angle-of-attack alignment. The bosses on each end of the upper airfoil 
(airfoil B) provided an attachment point to clamp the airfoil to the window (see fig. 8) at a fixed 0° angle 
of attack. The radius of the convex surface was 19.8 in. and the airfoil span was nominally 5.97 in., 
providing a 0.015-in. gap between each airfoil end and the test section window/side wall. The gap 
provided space for installing off-the-shelf felt pads with adhesive backing. One felt pad was placed 
upstream and one downstream of the trunnion on each end of the airfoil to prevent the airfoil from 
rubbing the test section window/side wall. The chord length was 5.94 in., representing a maximum 
thickness-to-chord ratio of 7.6 percent. The airfoils had knife leading and trailing edges. Although both 
leading and trailing edges had a nominal design half-angle of 8.7° for both the convex and wedge 
surfaces, the lower airfoil (airfoil A) and the upper airfoil (airfoil B) had inspected convex leading-edge 
half-angles of 8.7° and 8.9°, respectively. Seven midspan pressure ports having a diameter of 0.03 1 in. 
were provided at x/c locations of 0.12, 0.25, 0.375, 0.5, 0.628, 0.75, and 0.89. The pressure ports 
corresponding to these x/c locations were designated as ports A1 to A7 (leading edge to trailing edge) 
and ports B1 to B7 (leading edge to trailing edge) for airfoils A and B, respectively. 

Figure 12 shows that the airfoil is not symmetric about the chord line and that only one of the 
surfaces is convex. This asymmetry is possible because each cascade passage is aerodynamically isolated 
from all other passages in a cascade with an SLEL. Thus, these two monoconvex airfoils can simulate an 
infinite cascade of biconvex airfoils with an SLEL. 


3.3 Oscillation Mechanism 

The airfoil oscillation mechanism design (fig. 13) was inspired by the NASA Glenn Duct 
Laboratory’s 4- by 10-in. wind tunnel (ref. 23). The oscillation system consists of a barrel cam, a cam 
follower arm, and a cam follower button. Oscillatory motion about midchord was imparted to airfoil A by 
the rotating barrel cam (figs. 8 and 14) driven by a 5-hp electric motor (fig. 15). A circular cam follower 
button attached to the cam follower ami 3 in. away from the airfoil axis of rotation traveled in the three- 
cycle sinusoidal barrel cam groove. As the barrel cam rotated, the sinusoidal cam groove transmitted an 
oscillatory motion to the cam follower button, which in turn provided an oscillatory motion to the cam 
follower arm and thus to the airfoil. The cam follower arm and button were connected to the airfoil boss 
as shown in figure 16. The barrel cam was mounted on a shaft by using a split taper bushing. The mean 
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angle of attack about which airfoil A oscillated could be adjusted by the longitudinal position of the 
barrel cam on the shaft as indicated in figure 13(b). 

The cam and cam button were made of 440C stainless steel. The cam was heat treated to Rockwell 
62 hardness. The button was heat treated to Rockwell 55 hardness to ensure that the button would wear 
before the cam. The cam, follower arm, and button were mounted inside a Plexiglas lubrication enclosure 
(see fig. 8). Lubrication was provided by completely immersing the cam, follower arm, and button in 
85W-140 gear oil. 


4.0 Instrumentation 

The instrumentation used in the experimental facility falls into two groups. The first group, called 
steady-state data instrumentation, was used for acquiring the various wind tunnel calibration data, 
including the steady-state airfoil pressures and the wind tunnel operating conditions, for each wind tunnel 
run. The second group, called unsteady data instrumentation, was used to measure the unsteady pressures 
on the oscillating cascade. The following sections describe the steady and unsteady instrumentation. 


4.1 Steady-State Data Instrumentation 

A thermocouple located upstream of the last screen of the settling chamber was used to measure 
stagnation temperature (see fig. 4). A total pressure probe located downstream of the last screen in the 
convergent portion of the wind tunnel was used to measure stagnation or total pressure. The total 
pressure probe was connected in parallel to the Scanivalve pressure -measuring system and to a separate 
total pressure transducer. Static pressure was measured on the floor of the test section by using a static 
pressure port. The static port was located upstream of any possible oblique shock waves emanating from 
the airfoils. The static port was likewise connected in parallel to both the Scanivalve pressure transducer 
and a separate static pressure transducer. The separate transducers were used to monitor the tunnel total 
and static pressures throughout the wind tunnel run. 

The Scanivalve system consisted of two sample-and-hold modules connected to one Scanivalve 
pressure transducer with a range of 0 to 100 psig and a sensitivity of 0.2 mV/psi. The separate total and 
static pressure transducers had a range of 0 to 100 psig and a sensitivity of 0.01 V/psi. 

Bay signal conditioning units provided excitation voltage, bridge balancing, and amplification for the 
Scanivalve, the wind tunnel steady-state pressure transducers, and the thermocouple. The steady-state 
signals (except for the steady-state signals from the Kulite pressure transducers mounted in the airfoil) 
were passed through a low-pass filter set at 10 Hz. The steady-state pressure signals were converted to 
digital signals by using a ±1-V Datel model 256 analog-to-digital (A/D) convertor with 14-bit resolution. 


4.2 Unsteady Data Instrumentation 

Kulite model XCS-093-25-A pressure transducers having a range of 0 to 25 psia were installed in 
the airfoils. The transducers had a nominal sensitivity of 10 mV/psi at 5-Vdc excitation. Each transducer 
(with an outside diameter of 0.095 in.) was mounted in a stainless steel sleeve with an inside diameter of 
0.1 in. and an outside diameter of 0.13 in. The sleeve and transducer assemblies were in turn mounted in 
the 0.138-in.-diameter airfoil instrumentation orifices as shown in figure 17. Dow Coming 3140 room- 
temperature-vulcanizing (RTV) silicone coating was used to mount the transducer in the sleeve and the 
sleeve/transducer assembly in the airfoil instrumentation orifices. The sleeve provided a means for 
installing and removing the pressure transducers without damaging them. The sleeve could be gripped 
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with pliers, or a small wire could be inserted into a hole in the sleeve side wall to pull it out of the airfoil 
(design concepts from private communication with David Seidel of NASA Langley, February 1991). 

The transducers were mounted parallel to the airfoil span such that their pressure-sensing diaphragms 
lay in the plane of oscillation (chordwise direction), perpendicular to the airfoil axis of rotation. In this 
manner any apparent pressure response due to motion induced into the pressure diaphragms from the 
airfoil oscillation was minimized. 

The phase angle of the unsteady pressures with respect to the airfoil oscillatory motion were 
determined by using an Endevco model 7264-2000 accelerometer mounted on the follower arm. The 
accelerometer was encased in RTV to protect the sensor from the lubricating oil. At 10-Vdc excitation, 
the nominal output sensitivity was 0.2 mV/g. The accelerometer generated a sinusoidal signal as the 
airfoil oscillated. The phase angle was computed as the difference between the phase of the unsteady 
pressures and the phase of the airfoil motion. 

The excitation and bridge-balancing units used for the Kulite pressure transducers and the Endevco 
accelerometer had been manufactured by North American Aviation many years previously. Bay 5303 
signal-conditioning units were used to amplify and filter the analog signals. The excitation bridge- 
balancing and signal-conditioning units were all contained in a portable two-rack system as shown in 
figure 18. The analog signals were then converted to digital signals by using a Preston ±10-V A/D 
convertor with 15-bit resolution. 

Although the Kulites were absolute pressure transducers, they were used as differential transducers. 
Before each wind tunnel run the transducer bridges were balanced. The output voltage from the 
transducer (steady state) corresponded to the ambient atmospheric pressure. The unsteady pressures 
obtained were an offset from the atmospheric. Any damage to a transducer diaphragm would be apparent 
by inspecting the output voltage (steady state), with the correct voltage corresponding to atmospheric 
pressure. 

A crucial aspect of this type of experiment is to maintain a constant airfoil oscillation frequency. 
When obtaining unsteady pressures on a cascade with stagger by using the influence coefficient 
technique, two separate wind tunnel runs are required in which the airfoil oscillation frequencies need to 
be constant and the frequencies between the two runs must match. The airfoil oscillation frequency was 
monitored at the airfoil drive mechanism by using a 12-V digital tachometer connected to a magnetic 
probe. The magnetic probe was placed in close proximity to the drive pulley split taper bushing. The 
bushing contained a slot that provided a once-per-revolution signal as it passed by the probe. Airfoil 
oscillation frequency was determined from knowing the diameter ratio of the drive pulley to the cam 
pulley and taking into account that there were three airfoil oscillations to one revolution of the cam. This 
tachometer setup was used for determining airfoil oscillating frequencies and for matching the oscillating 
frequencies between wind tunnel runs for cascades with stagger. The airfoil frequencies used for 
calculating reduced frequencies were determined by performing a discrete Fourier transform (DFT) on 
the accelerometer signal. 

For unsteady testing the oscillation system was started before the wind tunnel. When steady-state 
airfoil oscillation had been achieved, the wind tunnel was started. Unsteady data were collected when the 
tunnel conditions indicated supersonic flow. 


5.0 Data Acquisition 

Steady-state pressures were acquired by using both the Scanivalve system and the Kulite pressure 
transducers. With the wind tunnel test section at Mach 1.9 the Scanivalve sample-and-hold modules were 
activated, trapping the pressures in the modules, where they were subsequently sampled in sequence by 
using the Scanivalve pressure transducer. All wind tunnel calibration data and the first series of steady- 
state cascade pressures were acquired in this manner. 
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As previously stated, the Kulites were also used to measure the steady pressures on the airfoils. The 
acquisition technique for steady pressure data using the Kulites was the same as for the acquisition of 
unsteady data described in the following text. 

The magnitude and phase of the unsteady pressures were acquired by using Kulite pressure 
transducers and an Endevco accelerometer as previously mentioned. Fifteen Preston sample-and-hold 
cards were used for simultaneous sampling of the 14 Kulite pressure transducer signals and the Endevco 
accelerometer signal. With the tunnel operating at Mach 1.9 and airfoil A oscillating at a near constant 
frequency, the electronic sample-and-hold cards were activated while retaining the various voltages 
output from the transducers until they were processed by the Preston A/D convertor. The data sample- 
and-hold process was accomplished at a rate of 10 kHz. After the unsteady pressure data were acquired, 
the wind tunnel total temperature and total and static pressures were recorded by using the Scanivalve 
system. The tunnel was then shut down. 

A wideband filter was used to attenuate signals above 25 kHz. Figure 19 shows a typical calibration 
curve from one of the amplifiers. The amplifiers provided a flat gain response from direct current to 
20 kHz, with a 3-dB rolloff at 25 kHz. According to the Nyquist frequency criteria, signals between 
5 kHz and 25 kHz would be “wrapped” around the DFT window. If any frequencies of these alias signals 
were an integer multiple of the unsteady pressure signals, erroneous pressure magnitudes would result. 
Upon inspection of a typical power spectrum of the unsteady pressure signals as shown in figure 20, one 
can intuitively see that there was insignificant signal content above the first harmonic. This result was 
typical. Thus, although aliasing did occur, it was not believed to be an overriding problem in our case. 


6.0 Wind Tunnel Calibration 

The flow characteristics of the wind tunnel were determined so that the unsteady aerodynamic data 
could be properly reduced. Mainly, wind tunnel calibration consisted of determining the test section 
Mach number and flow angle. These two quantities were determined prior to unsteady airfoil testing. 


6.1 Mach Number 

The test section Mach number was measured by operating the tunnel empty and measuring the static 
pressure on the roof and floor of the test section. Figure 21 is a plot of p/po for the floor and roof versus 
distance along the tunnel axis (starting after the throat and continuing through the test section). The test 
section p/po pressure ratio corresponded to that for Mach 1.97 flow. Thus, supersonic flow in the test 
section was confirmed. 


6.2 Flow Angle 

The next series of tests consisted of measuring the flow angle in the test section. A flow angle wedge 
(fig. 22) was installed in the airfoil A position in the tunnel test section. The flow angle wedge had one 
static pressure orifice at midspan on each side of its 18° (half-angle) wedge. The wedge was positioned at 
-1°, 0°, and 1° angles of attack. Comparing the upper and lower surface pressures determined that the 
flow was essentially parallel to the tunnel centerline at the test section. 

The presence of the airfoils in the wind tunnel would change the flow in the test section and could 
also introduce enough flow blockage to result in a subsonic leading-edge locus on the cascade. Thus, the 
flow field in the test section with the airfoils installed was investigated, in several steps, to check for any 
undesirable characteristics. First, the wind tunnel was operated with only airfoil A installed. Pressures on 
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the convex surface of the airfoil were measured at 1°, 0°, and -1° angles of attack. Figures 23(a) and (b) 
compare the experimental pressures obtained by using the Scanivalve system with pressures predicted by 
shock expansion theory for an airfoil at 0° and ±1° angles of attack, respectively. Experimental data at x/c 
of 0.625 were not included because an obstruction in the pressure tubing prevented proper pressure 
measurement. Experiment and theory correlated well. 

The wind tunnel was operated with both airfoils A and B installed. Steady-state pressures were 
measured on the convex surfaces of both airfoils at stagger angles of 0°, 10°, and 20°. Figure 24 compares 
experimental pressures (obtained by using the Scanivalve system) with Euler CFD code predictions 
(ref. 18) for a stagger angle of 0°. The correlation was very good. Figures 25 and 26 compare 
experimental pressures obtained by using the Scanivalve system and the Kulite pressure transducers with 
CFD prediction for stagger angles of 10° and 20°, respectively. The Kulites did not correlate as well with 
the theory at the trailing edge as did the pressures measured by the Scanivalve system. 

Looking at the correlation of theory with the Scanivalve data, it can be concluded that the cascade 
did indeed have a supersonic leading-edge locus. It can also be concluded that the free-stream Mach 
number was 1.97 and the flow behaved well. 


7.0 Unsteady Data Reduction 

This section describes the methods by which the unsteady experimental data were transformed into 
influence coefficients. The unsteady, analog, time -history voltage output from the transducers was 
converted into discrete time-history pressure and acceleration values by using an A/D convertor and the 
transducer sensitivities. Because the DFT was required to calculate the influence coefficients and the 

Yl . 

DFT routine (ref. 24) required 2 data points, the number of discrete data samples was limited to 

Yl 

quantities of 2 , where 1 < n < 15. For each transducer 32 768 (n = 15) data samples were acquired 
simultaneously, at a rate of 10 000 samples/sec, resulting in a data acquisition time of 3.2768 sec for each 
wind tunnel run. Consequently, the time increment between samples was 0.0001 sec. 

As previously stated, the 15 transducer time histories consisted of 32 768 data points each. The time 

Yl 

histories were then divided into blocks of data with 2 samples each, where n was either 11, 12, or 13, 
resulting in 2048, 4096, or 8192 points, respectively. The optimum choice of the order n was determined 
by calculating the unsteady pressure coefficients for all three values of n and then choosing the data with 
the smallest uncertainty based on a ±95 percent confidence interval obtained by using a “t” statistical 
distribution. 

One characteristic of a blowdown wind tunnel is that the total pressure and hence the static pressure 
in the test section decrease with time. Thus, the sinusoidal pressure -time history contained a 
nonoscillatory decreasing static pressure component. The data were consequently detrended to remove 
the nonoscillatory component such that the sinusoidal pressure -time history contained only unsteady 
pressure components. 

FORTRAN software was used to transform the 15 transducer time histories into 14 unsteady pressure 
magnitude and phase values. The code was provided by Dr. Dan Buffum of NASA Glenn and was 
modified for use at The Ohio State University. The flowchart in figure 27 presents, in a general and 
condensed manner, the process of transforming the transducer time histories into magnitude and phase. 
This process will be described next. 

Referring to figure 27, the accelerometer time history consisting of 32 768 data points was divided 
into 16 blocks of data (assuming n = 11) containing 2048 data points each. A Hanning window was 
applied to each block. The phase angle of the first harmonic of the airfoil oscillation signal for each 
block was determined by using a DFT routine and interpolating (ref. 25) to correct for leakage 
(see appendix B). The acceleration harmonics were converted into displacement harmonics before 
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interpolation. Thus, 16 phase angles of the airfoil oscillation, corresponding to the 16 data blocks, were 
determined. 

Next, the time history of the unsteady pressure for port A1 (the transducer nearest the leading edge of 
airfoil A) was divided into 16 data blocks containing 2048 data points each. A Hanning window was 
applied to each block. The magnitude and phase angle of the first harmonic of the unsteady pressure were 
determined by using a DFT and interpolating. However, of interest was the phase angle of the unsteady 
pressure with respect to airfoil oscillation, and not the phase angle of the unsteady pressure-time history 
itself. Therefore, the phase angle of interest was determined by subtracting the phase angle of the airfoil 
oscillation time history from the phase angle of the unsteady pressure-time history, for all 16 data blocks 
for port Al. 

An average magnitude and phase angle of the unsteady pressure -time history were determined by 
using all 16 data blocks. The magnitude and phase angle resulting from each data block were then 
compared against the average. Blocks of data with their magnitudes outside ±20 percent of the average 
were discarded. Blocks of data with their phase angles in a quadrant opposite the average were also 
discarded. A new average magnitude and phase angle were computed by using the remaining data blocks. 
This process was then repeated for the remaining pressure transducers on airfoils A and B. 

Thus, the magnitudes and phase angles (with respect to airfoil oscillation) of the first harmonic of the 
unsteady pressure for 14 transducers were acquired. From the magnitudes and phase angles the real and 
imaginary components were computed. Each unsteady pressure having a magnitude and phase angle is 
represented as the complex quantity P\ in equation ( 1 ). 

The influence coefficients were calculated by dividing the real and imaginary components of P\ by 
q.jSi i , where is given as 


(23) 

where y is the ratio of specific heats (y = 1.4), p yj is the free-stream static pressure, and M, x is the free- 
stream Mach number. 


7.1 Free-Stream Mach Number 

Assuming isentropic flow, the Mach number was determined by using the measured total pressure in 
the stagnation chamber and the measured static pressure on the test section floor. Although the total 
pressure and hence the static pressure decreased during each wind tunnel run, this decreased pressure did 
not affect the Mach number because it was fixed by the tunnel geometry (ratio of throat to test section 
area), given that the ratio of static to total pressure was sufficient to maintain supersonic flow. 


7.2 Free-Stream Static Pressure 

Because the free-stream static pressure poo decreased during a wind tunnel run, an average value was 
used in equation (23) to calculate the influence coefficients. The static pressure port used to determine 
Mach number was sampled only at the end of each wind tunnel run. Thus, using this port to provide the 
static pressure would have resulted in a lower than average static pressure in equation (23) and hence 
higher than average influence coefficients, which would not be desirable. 

Because airfoil B did not oscillate and port B1 was upstream of any shock reflections, the free-stream 
static pressure was determined by using shock expansion theory and the measured pressure from port B 1 . 
Figure 28 presents a flowchart describing this method. From the free-stream Mach number, the airfoil 
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surface slope at port B 1 , and the time history of the pressure p-Q from port B 1 , a time history of the free- 
stream static pressure was obtained. An average was determined from the detrended data and used to 
calculate qco- 


7.3 Amplitude of Oscillation 

The amplitude of oscillation was measured statically before and/or after the wind tunnel run by using 
a gunners sight (fig. 29). By clamping a steel plate to the top horizontal flat surface on the north airfoil 
boss and resting the gunners sight on the plate, the angle of the plate (and hence the airfoil chord line) 
with respect to level was measured. The angle was measured at the extreme pitch-up and pitch-down 
positions of the airfoil by contacting the cam follower button on the appropriate inside wall of the barrel 
cam groove. This measurement was made for all lobes of the barrel cam. Even though the gunners sight 
had an accuracy to within ±0.0056°, the measurement of the airfoil oscillation amplitude was not as 
accurate, in part because of minute deviations in the flatness of the plate. The method of clamping the 
plate to the airfoil also introduced errors. The wear of the cam follower button introduced play between it 
and the cam groove, which also made it difficult to measure the exact value of the oscillation amplitude. 
Despite uncertainties in the exact oscillation amplitude, an amplitude of 1.2° was about midway between 
the estimated error limits and therefore was used for data reduction. The pitch-up and pitch-down 
amplitudes were set as nearly equal as possible prior to the wind tunnel runs such that the airfoil 
oscillated approximately about 0° angle of attack. 


8.0 Unsteady Aerodynamic Results 

This section compares experimental unsteady pressures with those predicted by three unsteady 
aerodynamic algorithms for a cascade of airfoils with an SLEL oscillating in pitch about midchord. 

These algorithms were based on potential theory, piston/potential theory, and the Euler equations. 

Reference 9 presents a potential theory solution for the unsteady pressure on a cascade of oscillating 
airfoils with an SLEL. This approach idealizes the airfoils as flat plates. The streamwise location of 
shock wave/airfoil impingement is determined by using Mach angles. Nonlinear effects on unsteady 
pressure and shock wave/airfoil impingement location, due to the airfoil geometry, are not accounted for. 
Reference 8 presents an algorithm based on this theory, denoted herein as the Lane code. 

In an effort to maintain computational efficiency and include some effects of airfoil geometry, a 
hybrid algorithm was developed to include the nonlinear effects of thickness and camber in the isolated 
noninterference regions of the airfoils. These regions consisted of the portion of each airfoil upstream of 
the first shock wave/airfoil impingement location. Because of the SLEL each airfoil surface in this region 
was aerodynamically isolated from its neighbors. Piston theory was used in this region to account for the 
nonlinear effects of thickness and camber. Downstream of the first shock wave/airfoil impingement 
location the potential theory of reference 9 was used. Even though airfoil geometry was accounted for in 
the isolated noninterference regions of the airfoils, the shock wave/airfoil impingement locations were 
determined by assuming flat-plate airfoils in conjunction with the use of Mach angles. Reference 12 
presents the algorithm and theoretical development for this hybrid approach, denoted herein as the 
piston code. 

Reference 1 8 describes the computational fluid dynamic algorithm used to solve the Euler equations 
and provide unsteady pressures for comparison with experiment. This algorithm, referred to herein as the 
Euler code, solves the two-dimensional unsteady Euler equations by using a time-marching, flux- 
difference splitting scheme of reference 17 with blocked grids. The Euler code solves the flow equations 
for one or more El grids, which are allowed to deform with airfoil motion as in reference 16. The 
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algorithm is third-order accurate in space and second-order accurate in time. The Euler code accounts for 
the airfoil geometry while solving the inviscid flow equations. 

Because comparing experimental and predicted unsteady pressure coefficients for a variety of 
interblade phase angles would be a time-consuming task, the comparisons presented in this report are for 
interblade phase angles of 0° and 180° only. However, the reader may generate experimental unsteady 
pressure difference coefficients for any interblade phase angle by using the unsteady pressure influence 
coefficients presented in table 1 and assigning an interblade phase angle value to P in equation (15). 

Influence coefficient data covering interblade phase angles of 0° and 180° are presented for 0° and 
10° of stagger. The unsteady data (data set) are identified by the corresponding wind tunnel run numbers. 
Table 11 presents the wind tunnel operating conditions for the data sets. 


8.1 0° Stagger 

Figure 30 compares the experimental values of the real and imaginary components and the magnitude 
and phase angle of A C p for data set 68 with the values predicted by the Lane and piston codes. Figures 3 1 
and 32 present these comparisons and comparisons with the Euler code predictions for data sets 69 and 
77, respectively. Data sets 68, 69, and 77 correspond to reduced frequencies of 0.0562, 0.0929, and 
0.1457, respectively. Because of undetermined computational difficulties, the Euler code did not run 
at the low reduced frequency of 0.0562; thus, no comparison was made between its predictions and data 
set 68. 

Figures 30 to 32 show that the various components of the experimental values of A C p change 
abruptly near x/c = 0.5. This abrupt change resulted from the bow shock waves impinging on the 
downstream airfoil surfaces. The spacing of the seven pressure transducers made it impossible to 
determine the precise location of the shock wave reflection. For future reference the shock wave/airfoil 
impingement location from experiment is denoted as x s . 

8.1.1 Real and imaginary parts of A C p . — The Rc(AC /; ) changes from a positive quantity upstream 
of x s to a negative quantity downstream (figs. 30(a), 3 1(a), and 32(a)), but the lm(AC /; ) behaves in the 
opposite manner (figs. 30(b), 31(b), and 32(b)). Examining the Re(AC^) and lm(AC /; ) in these figures 
shows that the correlation between experiment and the various algorithms is not as favorable for Im(AC^) 
as for Re(ACp). The shock wave/airfoil impingement location predicted by the Lane and piston codes for 
lm(A Cp) is about 20 percent downstream of its actual location. For future reference, x s denotes the 
analytically predicted streamwise location of bow shock wave/airfoil impingement based on the Lane and 
piston codes. 

As shown in figures 31(a) and (b) and 32(a) and (b), the Lane code does not correlate with 
experiment as well as the Euler code for Re(AC /; ) and Im(ACp) upstream of x s . This result is probably 
due to airfoil surface curvature, which the Lane code does not account for. The piston code does a better 
job than the Lane code in this respect. The Euler code correlates better with experiment than either the 
Lane or piston codes downstream. 

8.1.2 Magnitude of AC p . — Figures 30(c), 31(c), and 32(c) compare experimental and predicted 
values of the magnitude of A C p , denoted as \AC p \ for data sets 68, 69, and 77, respectively. As shown in 
figures 3 1(c) and 32(c), the Euler code correlated better than the Lane and piston codes. The Euler code 
exhibited the same trends as the experiment, with the largest discrepancies occurring around midchord. 

8.1.3 Phase angle of A C p . — Figures 30(d), 31(d), and 32(d) compare experimental and predicted 
phase angles for data sets 68, 69, and 77, respectively. The Euler, Lane, and piston codes all correlated to 
the same degree with experiment upstream of x s , . Again, the Lane and piston codes predicted x s to be 
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20 percent farther downstream than experiment showed. The Lane and piston codes predicted the shock 
wave/airfoil impingement location to be farther downstream owing to their inherent assumptions as 
discussed previously. Downstream of x s the Lane and piston codes predicted phase angles that 
correlated reasonably well with the experimental values downstream of x s . The predicted phase angles 
from the Euler code correlated well with experiment downstream of x s . 

After inspection of the unsteady data for 0° stagger it was evident that the experimental and predicted 
values of the imaginary part of A C p did not correlate as well as the real part in the isolated airfoil regions. 
The source of this discrepancy is believed to be in part the difference between the experimental and 
analytical phase angles. As shown in figures 30(d), 31(d), and 32(d) and schematically in the complex 
plane in figure 33, the phase angles upstream of x s are small. Figure 33 shows that any error in the 
small phase angle a will affect the imaginary component to a greater extent than the real component. To 
further corroborate this idea, the phase angles calculated by using the Euler code were used in 
conjunction with the experimental unsteady pressure magnitudes of data set 69 to determine the real and 
imaginary components of the unsteady pressure difference A C p . The real and imaginary components of 
A C p correlated quite well with the Euler-code-calculated components of A C p in the isolated airfoil 
regions, as shown in figure 34. This correlation confirms that the measurement of the phase angle was 
one source of the discrepancy between experimental and analytical unsteady pressures. The correlation 
was also improved downstream of the shock wave/airfoil impingement location, but not as well as 
upstream. 

Several possible causes of phase angle error were possible, including airfoil twisting during 
oscillation, acoustic effects of the pressure orifice, excessive wear of the cam follower button, as well as 
possible accelerometer error. Although a malfunctioning accelerometer could have contributed to phase 
angle error, it was not investigated. However, the other three possible causes were investigated and are 
discussed next. Airfoil twisting during oscillation could change the actual phase angle between the 
accelerometer and the pressure transducers because the accelerometer was located on the airfoil arm and 
the pressure transducers were located at midspan of the airfoil. The phase angle error due to airfoil 
twisting during oscillation was determined to be negligible by performing a direct transient response on a 
finite element model of the airfoil using MSC/NASTRAN. The airfoil was modeled by using solid 
elements. An enforced sinusoidal (time domain) displacement was applied to the airfoil. Oscillation 
frequencies of 96.14 and 165.74 Hz were imposed. In both cases the phase lag due to airfoil twist was 
less than 0.4°. There was essentially no amplitude increase. 

The phase angle error due to pressure orifice geometry was investigated by Dr. Stephen Whitmore of 
NASA Ames/Dryden using a CFD code (ref. 26). The analysis estimated that the orifice geometry 
introduced a pressure phase lag of about 0.4° and a pressure rise of 0.05 dB. These values were 
considered negligible. 

Some phase angle error was very likely introduced by excessive wear of the cam follower button. As 
the cam follower button wore, the button would tend to slap inside the barrel cam groove instead of 
smoothly following the groove contour. Thus, the sinusoidal motion of the button and hence the airfoil, 
as well as the peak airfoil oscillation amplitude, were adversely affected. 

8.1.4 180° interblade phase angle . — The previous results for 0° stagger corresponded to an 
interblade phase angle of 0°. As mentioned before, the experimental influence coefficients may be used 
to generate the unsteady pressures for any interblade phase angle. As an example, the influence 
coefficients presented in table 1(a) are substituted into equation (15) with P = n, thereby generating the 
unsteady pressure difference for an interblade phase angle of 180°. Figures 35 to 37 show the real and 
imaginary parts of the unsteady pressure difference for runs 68, 69, and 77, with an interblade phase 
angle of 180°. Comparing these figures for P = 180° with figures 30(a) and (b), 3 1(a) and (b), and 32(a) 
and (b) for P = 0° shows that the unsteady pressure difference upstream of x s remained unaffected by 
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the interblade phase angle P , whereas downstream of x s the unsteady pressure difference changed both 
magnitude and sign. The correlation between experiment and the Lane code is lacking, although the Lane 
code exhibited the general trends. 


8.2 10° Stagger 

Figure 38 compares the experimental values for Re(ACp) and Im(ACp) as well as the magnitude and 
phase angle of A C p with the values predicted by the Euler, Lane, and piston codes for the cascade at 10° 
stagger. Data sets 84 and 86 were used to generate the influence coefficients. Data set 84 represents 
influence coefficients C® and C~ 1 for a stagger angle of 10°. Data set 86 represents influence 
coefficients C® and C l p for a stagger angle of -10°. Together these two data sets were used to produce 
the graphs shown in figure 38. 

Table 11 presents the various tunnel flow data for runs 84 and 86. The difference in reduced 
frequency between data sets 84 and 86 was less than 5.6 percent. The average ffee-stream static pressures 
and Mach numbers were within 1 percent and 0.56 percent of each other, respectively. 

8.2.1 Real and imaginary parts of A C p . — As shown in figure 38(a) the Euler code correlated better 
with the Re(ACp) than either the Lane or piston codes. Again, the Lane and piston codes placed the shock 
wave/airfoil impingement locations farther downstream than did the experiment for the aforementioned 
reasons. Figure 38(b) shows that the experimental values of the Im(AC),) did not correlate well with any 
of the codes. The lm(AC 0 ) as predicted by the Lane and piston codes would correlate better between x/c 
values of 0.4 and 0.78 if the analytical shock wave/airfoil impingement locations were shifted upstream 
by 20 percent. 

8.2.2 Magnitude of A C p . — Figure 38(c) compares the experimental and calculated values of \AC p \. 
Again, the inherent assumptions of the Lane and piston theories placed the analytical pressure jumps due 
to shock wave/airfoil impingement approximately 20 percent farther downstream than did the 
experiment. 

8.2.3 Phase angle of A C p . — As shown in figure 38(d) the phase angle as predicted by the Euler code 
correlated better with experiment than did the Lane or piston codes. The Lane, piston, and Euler codes 
were in close agreement upstream of 0.4 x/c. 

8.2.4 180° interblade phase angle. — The previous results for the 10° stagger cascade corresponded 
to an interblade phase angle of 0°. The unsteady pressure difference for an interblade phase angle of 180° 
was generated by substituting the influence coefficients from table 1(b) into equation (15) with P = n. 
Figure 39 shows the real and imaginary parts of the unsteady pressure difference for runs 84 and 86 with 
an interblade phase angle of 180°. Comparing these figures for P = 180° with figures 38(a) and (b) for 
P = 0° shows that the unsteady pressure difference upstream of the first shock wave/airfoil impingement 
location (x/c < 0.375) remained unaffected by the interblade phase angle. Correlation between 
experiment and the Lane code appeared to be lacking. 


9.0 Conclusions 

An experimental influence coefficient technique was used to obtain unsteady aerodynamic influence 
coefficients and unsteady pressures for a cascade of airfoils oscillating in pitch about midchord. Stagger 
angles of 0° and 10° were investigated for a cascade with a gap-to-chord ratio of 0.417 operating at an 
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axial Mach number of 1.9, resulting in a supersonic leading-edge locus. Reduced frequencies ranged 
from 0.056 to 0.2. Influence coefficients were obtained from which the unsteady pressures for any 
interblade phase angle could be determined. The unsteady pressures were compared with those predicted 
by several algorithms for interblade phase angles of 0° and 180°. 

This experimental program has laid the foundation for future unsteady aerodynamic cascade testing 
in the supersonic leading-edge -locus flow regime. The correlation of experimental and analytical results 
was sufficient to conclude that the two-airfoil technique is a good approach for simulating an oscillating 
cascade of symmetric airfoils with a supersonic leading-edge locus and obtaining experimental influence 
coefficients, ft is hoped that this effort will be of value to future experimental unsteady aerodynamic 
testing endeavors. 

The following conclusions pertain to this specific experimental effort: 

1 . Piston theory does reflect the general effects of airfoil geometry on unsteady pressure in the 
isolated airfoil regions of this cascade. 

2. The Lane, piston, and Euler codes all predict similar phase angles between the airfoil motion and 
the unsteady pressure. Thus, it is evident that the airfoil thickness distribution does not have a 
large effect on the unsteady-pressure phase angle. 

3. In general, the Euler code correlates better with experiment than do either the Lane or piston 
codes for the airfoils considered in this study. 

4. One major finding of this work is that the Lane and piston codes do not accurately locate the 
shock wave/airfoil impingement locations partly because of flat-plate airfoil assumptions. Even 
though the piston code accounts for airfoil geometry in the isolated airfoil regions, the shock 
wave/airfoil impingement locations are still based on flat-plate geometry. The unsteady pressures 
predicted by the Lane and piston codes could be greatly improved by repositioning the shock 
wave/airfoil impingement locations by using a steady-flow theory (e.g., the method of 
characteristics or a steady CFD flow solver) that accounts for airfoil geometry. 

5. Phase angle errors can be greatly reduced in future endeavors by improving the airfoil/ 
oscillating mechanism. These improvements can be incorporated with little difficulty. In 
particular, the airfoils could be made from a lighter material, such as titanium, along with smaller 
trunnions, reducing the airfoil/oscillating mechanism inertia characteristics and thereby reducing 
cam follower button wear. The cam follower button diameter could also be made larger to 
provide better wear performance. Minimizing cam follower button wear would provide more 
accurate sinusoidal airfoil oscillations, resulting in more accurate phase angles of the unsteady 
aerodynamic influence coefficients and thus providing more accurate unsteady pressures. 
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Appendix A 
Symbols 


A 

A g 

A h 

A r 

b 


C 


Pi 


c 


c 


Pu 

0 


Pi 


c 1 


Pi 


c~ 


Pi 


c u 


Pu 


c 1 


Pu 


C _1 

Pu 

A C p 
c 
Cl 
C m 


f 

fg 

fu,u+l 

g 

Im(. . .) 

k 

l 

M Uj u + 1 


Mqo 

N 


calculated estimate of amplitude 
given amplitude of cosine wave 

maximum amplitude of the weighting function for Hanning window 
maximum amplitude of the weighting function for rectangular window 
half chord, c/2 

pressure coefficient for lower surface of reference airfoil (n = 0) as defined in eq. (2) 
pressure coefficient for upper surface of reference airfoil (n = 0) as defined in eq. (2) 

influence coefficient for unsteady pressure on lower surface of reference airfoil (n = 0) 
due to oscillation of reference airfoil (see fig. 1) 

influence coefficient for unsteady pressure on lower surface of reference airfoil (n = 0) 
due to oscillation of airfoil n= 1 (see fig. 1) 

influence coefficient for unsteady pressure on lower surface of reference airfoil (n = 0) 
due to oscillation of airfoil n = — 1 (see fig. 1) 

influence coefficient for unsteady pressure on upper surface of reference airfoil (n = 0) 
due to oscillation of reference airfoil (see fig. 1) 

influence coefficient for unsteady pressure on upper surface of reference airfoil (n = 0) 
due to oscillation of airfoil n = 1 (see fig. 1) 

influence coefficient for unsteady pressure on upper surface of reference airfoil (n = 0) 
due to oscillation of airfoil n = — 1 (see fig. 1) 

pressure difference coefficient as defined in eq. (14) 
airfoil chord 

two-dimensional lift coefficient as defined in eq. (12) 

two-dimensional moment coefficient as defined in eq. (13) 

sampling frequency 

calculated estimate of frequency 

given frequency of cosine wave 

frequency component from discrete Fourier transform corresponding to bin u or u + 1 
gravitational constant 
imaginary part of (. . .) 
reduced frequency, co b/V 

subscript or superscript denoting lower surface of reference airfoil (n = 0) 

amplitude of frequency component from discrete Fourier transform corresponding to 
bin u or u + 1 

free-stream Mach number 

(number of airfoils - l)/2; number of data samples 
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n 

P\ 

P> Psi 
Pb 

Poo 

Po, P U 

Q GO 

7? 

Re(...) 

s 

t 

u 

V 

X 

X M 


y 

a\ 

P 

Y 

8 

0 

% 

®u,u + 1 

v 
a 
C >a 

Gp 

co 


airfoil index; also order of discrete Fourier transform 
complex-valued first harmonic of unsteady pressure signal 

static pressure 

static pressure at airfoil port B 1 
tree-stream static pressure 

stagnation or total pressure, i = 1,2 
tree-stream dynamic pressure 

ratio of two largest magnitudes of frequency spectrum 
real part of (...) 

blade spacing as shown in fig. 1 
time 

subscript or superscript denoting upper surface of reference airfoil (n = 0); discrete Fourier 
transform bin index 
tree-stream velocity 
streamwise coordinate 

streamwise location of shock wave/airfoil impingement from neighboring airfoil, as predicted 
by algorithms 

shock wave/airfoil impingement location on lower surface of reference airfoil 
shock wave/airfoil impingement location on upper surface of reference airfoil 
shock wave/airfoil impingement location from analysis 

shock wave/airfoil impingement location from experiment 

transverse coordinate 

amplitude of reference airfoil oscillation 

interblade phase angle; frequency change between consecutive bins 
ratio of specific heats (1.4 in this work) 
surface slope 

stagger angle; calculated estimate of phase angle; oblique shock wave angle 
given phase angle of cosine wave 

phase angle of frequency component from discrete Fourier transform corresponding 

to bin u or u + 1 

ratio defined in eq. (B2) 

phase angle 

phase angle of airfoil acceleration-time history 

phase angle of airfoil oscillation (displacement)-time history 

phase angle of unsteady pressure -time history 

phase angle of unsteady pressure with respect to airfoil oscillation 

circular frequency 
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Appendix B 
Leakage Correction 

This appendix describes the phenomenon of leakage by example, which it is hoped will give the 
reader a basic understanding of leakage and how it is corrected for in data reduction procedures. The 
equations used for leakage correction in this appendix were obtained from the notes of Dr. Don Braun of 
NASA Glenn. Appendix A of reference 21 and reference 25 discuss the derivation of these equations, 
although different variables are used. 

C onsider a continuous cosine wave having an amplitude of A g = fra frequency off g = 1 000 Hz, and 
a phase angle of 0 g = 0°. Suppose this wave is sampled at a rate of F s = 4000 samples/sec, capturing 
N= 4096 discrete data samples. Furthermore, suppose a rectangular window (each data point multiplied 
by a factor of A r = 1) is applied to the sampled data starting at t = 0. The end of the rectangular window 
would be at N/F s = 1.024 sec, corresponding to the total time for data acquisition. Exactly f g N/F s = 1024 
cycles of the cosine wave would be within the rectangular window, represented by the 4096 discrete data 
samples. 

Performing a discrete Fourier transform on the 4096 data samples within the rectangular window 
would result in the frequency spectrum shown in figure 40(a). There all components in the frequency 
domain have an amplitude of 0 except for the frequency components at +1000 Hz. The amplitudes of the 
frequency components at ±1000 Hz are both 0.5, each corresponding to one-half of the amplitude of the 
continuous cosine wave. 

It is sufficient to look only at the positive frequency spectrum and multiply the amplitude of the 
frequency component at 1000 Hz by 2 to get the correct amplitude of the cosine wave. This result is 
possible only because there was an integer number of cosine wave cycles in the rectangular window. 
Consequently, a frequency component at 1000 Hz in the frequency spectrum corresponded to the exact 
oscillation frequency (1000 Hz) of the cosine wave. Thus, multiplying the amplitude of the 1000-Hz 
component by 2 produced the correct amplitude. Moreover, inspecting the actual real and imaginary 
components of the amplitude of the 1000-Hz frequency component would show that the phase angle is 
0°. However, rarely is there an integer number of cycles of a signal in a data window. When the data 
window does not contain an integer number of cycles of a signal, “leakage” will result. 

C onsider a continuous cosine wave having an amplitude of A g = fra frequency of f g = 1 000 Hz, and 
a phase angle of () g = 0°, just as before. However, suppose this wave is sampled at F s = 3900 
samples/sec, capturing N = 4096 discrete data samples. Furthermore, suppose a rectangular window is 
also applied to the sampled data starting at t = 0 sec. The end of the window would be at N/F s = 1.050256 
sec (to seven significant digits), corresponding to the total time for data acquisition. There would be 
approximately f g N/F s = 1050.256 cycles of the continuous cosine wave within the window, represented 
by the 4096 discrete data samples. 

Performing a discrete Fourier transform on the 4096 data samples within the rectangular window 
would result in the positive frequency spectrum shown in figure 40(b) where each DFT frequency bin has 
a width of P = F s /N= 0.9521484 Hz. A closer look at the positive frequency spectrum (fig. 40(c)) reveals 
no frequency component at the cosine wave frequency of 1 000 Hz. The same situation occurs in the 
negative frequency spectrum, both the result of not having an integer number of cosine wave cycles in 
the data window. The frequency components containing the largest amplitudes are located in bin number 
u = 1050 at frequency f, = «[3 = 999.7559 Hz and in bin u + 1 = 1051 at f u +\ = (u + 1)|3 = 1000.708 Hz. 
Taking power from the 1000-Hz frequency component that existed in the continuous cosine wave and 
distributing this power to frequency components in the frequency spectrum that did not exist in the 
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continuous wave, such as at 999.7559 Hz and 1000.708 Hz (the predominant components among many 
others), is known as leakage. 

Referring to figure 40(c), the amplitudes for the frequency components at f, = 999.7559 Hz and 
f u +\ = 1000.708 Hz ar e M u = 0.4477416 andM J(+ i = 0.1542751, respectively. The corresponding phase 
angles are 9„ = 46.14598° (0.8053993 rad) and Q u +\ = -133.8233° (-2.335657 rad), respectively. It is 
obvious that the correct frequency, magnitude, and phase angle of the continuous cosine wave do not 
correspond to any frequency components in the spectrum but lie between the frequency components with 
the largest amplitudes. The correct amplitude would be as represented by the dashed vertical line at 
1000 Hz. To approximate the correct amplitude and phase angle of the original continuous cosine wave, 
a method of interpolation is used. 

As previously stated, a rectangular window was applied to the data. Equations (Bl) to (B9) are used 
to correct for leakage when a rectangular window has been applied to the data. The following discussion 
presents the interpolation calculations for the 1050.256 cycles of data in the rectangular window. The 
first step in correcting for leakage is to compute the ratio of the two largest magnitudes of the frequency 
spectrum, denoted as M u and M u + \ (fig. 40(c)). The ratio yields 


R = 


M„ 


M 


U~\~ 1 


0.4477416 

0.1542751 


2.902229 


(Bl) 


Next the following ratio is computed: 


v»— ^ = 0.2562638 (B2) 

1 + R 

Dividing the sampling frequency by the number of discrete data samples gives the frequency change 
between consecutive DFT bins as 


F 3900 

B = — = = 0.9521484 Hz (B3) 

N 4096 

Because the frequency component at /,, has a higher energy content than that at f u +\, corresponding to 0 < 
v < 1/2, equations (B4) to (B6) are used to estimate the frequency, amplitude, and phase angle, 
respectively, of the cosine wave. If the frequency component at f u+ \ had a higher energy content than that 
at f u , corresponding to 1/2 < v < 1, equations (B7) to (B9) would be used instead. The estimate of the 
frequency is then given as 


f = fu+ V P 


(B4) 


The estimate of the amplitude is given as 

2 A sin 


( 7tV^ 


A : 


N 


v J 


A r sin(7tv) 


M„ 


2nv 
A r sin(Trv) 


-M„ 


(B5) 


where nv is in radians. The second part of equation (B5) is valid for large N, where A r is the weighting 
function for the rectangular window and is equal to 1 for the DFT program used in this example. 
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Finally, the estimate of the phase angle is given as 


e*e 


u 


n(N - l)v 
N 


(B6) 


where 2n should be added to that value of 0 if it is less than -n to make the final value of 0 in the 
standard range between -n and n. Upon substituting f„ M u , and 0„, which are known from the DFT, and 
v and P from equations (B2) and (B3) into equations (B4) to (B6), the following estimates of frequency, 
amplitude, and phase angle for the cosine wave are obtained: 

/ = 999.9999 Hz 
A = 1.000068 

0 = 0.0005194 rad = 0.02975831° 

As is shown, good approximations to the actual frequency, amplitude, and phase angle result. 

As previously stated, when 1/2 < v < 1 and a rectangular window is applied to the data, the following 
equations should be used instead of equations (B4) to (B6). 


f = f u+ 1 -(i-v)P 


(B7) 


2 A sin 


At 


7t(l-v) 

N 


J r sin[Tt(l-v)] 


—M, 


u + 1 


27t(l-v) 


A r sin[Tt(l-v)] 


-M, 


W+l 


(B8) 


and 


0«0 


W+l 


, n (N - 1)(1 - v) 
N 


(B9) 


where 2n should be subtracted from that value of 0 if it is greater than n. 

Because of its desirable characteristics, a Hanning window was used in the data reduction for this 
report. Presented here are the equations that should be used to correct for leakage when a Hanning 
window is applied to the discrete data. Equations (BIO) to (B12) should be used to approximate the 
frequency, amplitude, and phase angle of a waveform when 0 < v < 1/2: 

/ = /„+v P (B 10 ) 


2N\ 


1- 


A : 


TCV 

tan — 

u 

tanf— 

u 


. ( 2nvP 
ism 

l N J 


A h sin(rcv) 




A h sin(Tcv) 


(B 11) 


and 
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(B12) 


0 ~ 0 M - nv 

where 2rc should be added to the value of 0 if it is less than -tc. 

The following equations should be used to correct for leakage when a Hanning window is applied to 
the discrete data and 1/2 < v < 1 : 


f = fu+\ — (! - v)P 


( B 13) 


tan 


2N< 


1- 


A : 


7t(l - V) 

N 


tanf— 

u. 


. ( 2n(l-v)^ 

Ism — - - 

l N 


A h sin[7r(l-v)] 


47rv(2-v)(l-v) 

' ■ x ’ m u+ i 
A h sm^rcv) 


(B14) 


and 


0 * 0 1<+1 + tt( 1 - v) 


( B 15) 


where 2n should be subtracted from that value of 0 if it is greater than tc. 
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TABLE I.— UNSTEADY PRESSURE INFLUENCE COEFFICIENTS 
(a) Stagger angle 9 = 0° 


(al) Upper surface 


Ratio of 
streamwise 
coordinate to 
airfoil chord, 
x/c 

Run 

! 68 

i 69 

: 77 i 

j Reduced frequency, k \ 

j 0.0562 

I 0.0929 

| 0.1457 1 

Influence coefficients for unsteady 
pressure on upper surface due to 

oscillation of reference airfoil, 

Pu 

0.120 

-0.8154 

+0.1451/ 

-0.8528 

+0.1522/ 

-0.9122 

+0.2491/ 

0.250 

-0.7505 

+0.1201/ 

-0.7650 

+0.1174/ 

-0.8196 

+0.2032/ 

0.374 

-0.6887 

+0.1181/ 

-0.7067 

+0.1205/ 

-0.7383 

+0.2022/ 

0.500 

-0.7211 

+0.0870/ 

-0.7569 

+0.06557/ 

-0.7577 

+0.1126/ 

0.627 

-0.7963 

+0.1099/ 

-0.8350 

+0.09825/ 

-0.8837 

+0.1645/ 

0.750 

-0.6437 

+0.0919/ 

-0.6587 

+0.08151/ 

-0.6937 

+0.1350/ 

0.890 

-0.4829 

+0.05467/ 

-0.4883 

+0.03536/ 

-0.4893 

+0.03686/ 


Influence coefficients for unsteady 
pressure on upper surface due to 

oscillation of airfoil n = 1, c\ 

0.120 




_ 



0.250 




_ 



0.374 




_ 



0.500 

1.204 

-0.3005/ 

1.242 

-0.3774/ 

1.274 

-0.6315/ 

0.627 

1.404 

-0.3692/ 

1.402 

-0.4570/ 

1.536 

-0.9019/ 

0.750 

1.067 

-0.2937/ 

1.078 

-0.3751/ 

1.037 

-0.5857/ 

0.890 

0.6880 

-0.1744/ 

0.6819 

-0.2109/ 

0.6436 

-0.3052/ 


(a2) Lower surface 


Ratio of 
streamwise 
coordinate to 
airfoil chord, 

xlc 

Run 

; 68 

69 

! 77 | 

! Reduced frequency, k 

j 0.0562 

0.0929 

| 0.1457 | 

Influence coefficients for unsteady 
pressure on lower surface due to 

oscillation of reference airfoil, 

Pi 

0.120 

0.8154 

-0.1451/ 

0.8528 

-0.1522/ 

0.9122 

-0.2491/ 

0.250 

0.7505 

-0.1201/ 

0.7650 

-0.1174/ 

0.8196 

-0.2032/ 

0.374 

0.6887 

-0.1181/ 

0.7067 

-0.1205/ 

0.7383 

-0.2022/ 

0.500 

0.7211 

-0.0870/ 

0.7569 

-0.06557/ 

0.7577 

-0.1126/ 

0.627 

0.7963 

-0.1099/ 

0.8350 

-0.09825/ 

0.8837 

-0.1645/ 

0.750 

0.6437 

-0.0919/ 

0.6587 

-0.08151/ 

0.6937 

-0.1350/ 

0.890 

0.4829 

-0.05467/ 

0.4883 

-0.03536/ 

0.4893 

-0.03686/ 


Influence coefficients for unsteady 
pressure on lower surface due to 

oscillation of airfoil n = -1, C ,} 

Pi 

0.120 




_ 



0.250 




_ 



0.374 




_ 



0.500 

-1.204 

+0.3005/ 

-1.242 

+0.3774/ 

-1.274 

+0.6315/ 

0.627 

-1.404 

+0.3692/ 

-1.402 

+0.4570/ 

-1.536 

+0.9019/ 

0.750 

-1.067 

+0.2937/ 

-1.078 

+0.3751/ 

-1.037 

+0.5857/ 

0.890 

-0.6880 

+0.1744/ 

-0.6819 

+0.2109/ 

-0.6436 

+0.3052/ 
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TABLE I.— Concluded 

(b) Stagger angle 8 = 10°; runs 84 and 86; reduced frequency k = 0.2; Mach number M « 1.95 


Ratio of streamwise 
coordinate to 
airfoil chord, 

xtc 

Influence 
coefficients for 
unsteady pressure on 
lower surface due to 
oscillation of 
reference airfoil, 

C° 

Pi 

Influence 
coefficients for 
unsteady pressure on 
lower surface due to 
oscillation of airfoil 
n = -1, 

c _1 

P, 

Influence 
coefficients for 
unsteady pressure on 
upper surface due to 
oscillation of 
reference airfoil, 

C° 

Pu 

Influence 
coefficients for 
unsteady pressure on 
upper surface due to 
oscillation of airfoil 
n = 1, 

c 1 

Pu 

0.120 

0.9808 


-0.9686 





-.45781 


+.52271 


0.250 

0.8708 


-0.8645 





-.35921 


+.41301 


0.374 

0.7964 

-0.3817 

-0.7833 



-.35851 

+.70191 

+.40651 


0.500 


-1.201 

-0.6718 




+1.1151 

+.27351 


0.627 

0.9671 

-1.131 

-0.9913 



-.38181 

+1.2721 

+.29721 


0.750 


-0.7515 

-0.7418 




+.76541 

+.27171 


0.890 


-0.5506 

-1.814 

1.732 



+.45301 

+3.1501 

-1.5761 


TABLE II. — WIND TUNNEL OPERATING CONDITIONS 


Data set 

Stagger angle, 
0, 

deg 

Mach number, 
M 

Average free-stream 
static pressure, 
psia 

Reduced 

frequency, 

k 

68 

0 

am. jg -■ 

7.2988 


69 

0 


6.7356 

iflKRHH 

77 

0 


5.7019 

.146 

84 

10 


6.2651 

.200 

86 

-10 


6.2266 

.190 
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Figure 1 . — Cascade of airfoils. 



Figure 2. — Cascade of airfoils with supersonic leading-edge 
locus (SLEL). 
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Figure 3. — Comparison of three-airfoil and two-airfoil techniques. 
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Figure 3. — Concluded. Comparison of three-airfoil and two-airfoil techniques. 
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Figure 5. — Air supply for OCF. 
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Figure 6. — OCF test section with north wall and north nozzle block removed. 



Figure 7. — Nozzle block and airfoil window. 



Figure 8. — Airfoil oscillating mechanism. Click this caption in the online PDF 
version to view a movie of a typical wind tunnel run with an oscillating 
mechanism. 
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Figure 9. — Monoconvex airfoil. 
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Figure 10. — Cascade airfoil with instrumentation cover removed, showing seven unsteady pressure transducers installed. 



Figure 11 . — Outer surface of airfoil, showing instrumentation cover installed. 



Figure 12. — Supersonic cascade. 
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Sintered bronze 
bushings 




Position of cam longitudinally on shaft 
will provide adjustment of airfoil mean 
angle of attack about which it oscillates 

Figure 13. — Oscillating mechanism design, (a) View looking toward stagnation 
chamber, (b) Top view. 
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Barrel cam 


Split taper bushing 


Figure 14. — Barrel cam with three-cycle-per-revolution sinusoidal groove. 



Figure 15. — Oscillating mechanism assembly. 
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Figure 16. — Cascade airfoil, showing assembly of follower arm and cam 
follower button. 


NASA/TM- 


- 2004-211820 


40 




I 



NAS A/TM— 2004-2 11820 


41 


Gain, dB 



Figure 18. — Conditioning units for Kulite pressure 
transducers and Endevco accelerometer. (Voltage 
excitation and bridge-balancing units are located 
in left rack. Signal amplification and filtering units 
are located in right rack.) 


50 — 

o— 


400— O O-O-O- -0-0-0 -O' 


30 — 


Gain 

10 

100 

200 

O 0-00 00 - 0-0 


o 


20 < 

10 


'(3 O -O 0-0— 0-0-0 “0~0< 


o 


o 1 - 

10 ° 


I I ililil I I I Mill I i I i 1 1 ill I I I I 

101 102 103 104 4X104 

Frequency, Hz 



0 1000 2000 3000 4000 5000 

Frequency, Hz 


Figure 1 9. — Typical amplifier calibration curve. 


Figure 20. — Typical unsteady pressure amplitude 
spectrum. 
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Figure 21 . — Static pressure distribution for run 16. 
Free-stream Mach number M x = 1 .9702. 



Figure 22. — Flow angle wedge. 
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Figure 23. — Steady-state surface pressures for 
airfoil A at Mach 1 .97. (a) 0° angle of attack, 
(b) +1 ° angle of attack. 
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Figure 24. — Steady-state surface pressures for 
cascade at 0° stagger angle and Mach 1 .9587. 
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O Upper surfaced _ 

□ Lower surface) Ex P eriment 



airfoil chord, x/c 


A Upper surfaced _ 

O Lower surface) Ex P eriment 



airfoil chord, x/c 

Figure 25. — Steady-state surface pressures for 
cascade at 10° stagger angle (free-stream Mach 
number = 1 .97; ratio of blade spacing to air- 
foil chord s/c = 0.417. (a) With Scanivalve system 
(b) With Kulite transducers. 


O Upper surfaced 
□ Lower surface) Ex P er ' ment 



airfoil chord, x/c 


A Upper surface'! 
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Figure 26. — Steady-state surface pressures for 
cascade at 20° stagger angle (free-stream Mach 
number M^ = 1 .96; ratio of blade spacing to air- 
foil chord s/c = 0.417. (a) With Scanivalve system 
(b) With Kulite transducers. 
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Figure 27. — Process of transforming transducer time histories into magnitude and phase (e.g., 32 768 data 
points assumed here). 
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Figure 28. — Process for calculating free-stream static pressure. 
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Figure 29. — Gunners sight. 
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Figure 30. — Unsteady pressure difference for run 68 with interblade phase angle p = 0° (reduced 
frequency k = 0.0562; stagger angle 0 = 0°; free-stream Mach number = 1 .9446; ratio of blade 
spacing to airfoil chord s/c = 0.417). (a) Real part, (b) Imaginary part, (c) Magnitude, (d) Phase angle. 
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Figure 31 . — Unsteady pressure difference for run 69 with interblade phase angle p = 0° (reduced 
frequency k = 0.09288; stagger angle 0 = 0°; free-stream Mach number M ^ = 1 .9573; ratio of blade 
spacing to airfoil chord s/c = 0.417). (a) Real part, (b) Imaginary part, (c) Magnitude, (d) Phase angle. 
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Figure 32. — Unsteady pressure difference for run 77 with interblade phase angle p = 0° (reduced 
frequency k = 0.14577; stagger angle 0 = 0°; free-stream Mach number M ^ = 1 .9617; ratio of blade 
spacing to airfoil chord s/c = 0.417). (a) Real part, (b) Imaginary part, (c) Magnitude, (d) Phase angle. 


NAS A/TM— 2004-2 1 1 820 


51 



Im (A C p ) 


lACpI cos (t = Re (ACp) 
lACpI sin o- = lACpI o = Im (ACp) 


Figure 33. — Effects of phase angle error. (Assume 
that o correc t = 10°. If (^measured is 5° (50 percent 
off), then Im (AC p ) will be 50 percent off. However, 
Re (ACp) will only be 1 .1 5 percent off, since 
cos (10°) = 0.98481 and cos (5°) = 0.99619.) 
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Figure 34. — Unsteady pressure difference for run 69 
using phase angles from Euler code (reduced fre- 
quency k = 0.09288; interblade phase angle p = 0°; 
stagger angle 0 = 0°; free-stream Mach number 
= 1 .9573; ratio of blade spacing to airfoil chord 
s/c = 0.41 7). (a) Real part, (b) Imaginary part. 

(c) Magnitude. 
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Figure 35. — Unsteady pressure difference for run 68 
with interblade phase angle 3 = 180° (reduced 
frequency k = 0.0562; stagger angle 0 = 0°; free- 
stream Mach number M x = 1 .9446; ratio of blade 
spacing to airfoil chord s/c = 0.417). (a) Real part, 
(b) Imaginary part. 


Figure 36. — Unsteady pressure difference for run 69 
with interblade phase angle 3 = 180° (reduced 
frequency k = 0.09288; stagger angle 0 = 0°; free- 
stream Mach number M x = 1 .9573; ratio of blade 
spacing to airfoil chord s/c = 0.417). (a) Real part, 
(b) Imaginary part. 
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Figure 37. — Unsteady pressure difference for run 77 
with interblade phase angle p = 180° (reduced 
frequency k = 0.14577; stagger angle 6 = 0°; free- 
stream Mach number M x = 1.9617; ratio of blade 
spacing to airfoil chord s/c = 0.417). (a) Real part, 
(b) Imaginary part. 
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Figure 38. — Unsteady pressure difference for runs 84 and 86 with interblade phase angle p = 0° and stagger 
angle 0 = 10° (reduced frequency k = 0.2; free-stream Mach number M ^ = 1 .9532; ratio of blade spacing 
to airfoil chord s/c = 0.417). (a) Real part, (b) Imaginary part, (c) Magnitude, (d) Phase angle. 
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Figure 39. — Unsteady pressure difference for runs 
84 and 86 with interblade phase angle f$ = 180° 
and stagger angle 0 = 10° (reduced frequency 
k = 0.2; free-stream Mach number M x = 1 .9532; 
ratio of blade spacing to airfoil chord s/c = 0.417), 
(a) Real part, (b) Imaginary part. 


Figure 40. — Frequency spectrum of continuous 
cosine wave of amplitude 1 and frequency of 
1 000 Hz using rectangular data window. 

(a) Integer number of cycles in data window. 

(b) Noninteger number of cycles in data 
window, (c) Detailed view of (b). 
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